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1. Introduction 

The fuzzy approximation scheme [[T]] consists in approximating the algebra of functions on a 
manifold with a finite dimensional matrix algebra instead of discretising the underlying space as a 
lattice approximation does. 

Here we report our results for a hermitian scalar field on the fuzzy sphere. We find the col- 
lapsed phase diagram and in particular we calculate the uniform ordered/non-uniform ordered line 
that was absent in |^. 

The cuiTcnt study could be relatively easily repeated for a hermitian scalar field on other fuzzy 
spaces. The simplest extension would be to fuzzy CP^. Some variants of the scheme can be applied 
to fuzzy versions of and S'^^. The study reveals that the non-uniform disordered phase lines 
should correspond to a pure matrix model transition. 

As an approximation scheme, this "fuzzification" is well suited to numerical simulations of 
field theories [Q]. As a test run, the first fuzzy approximation to be investigated should be the 
simplest one, that of the two dimensional sphere CP^ = S^. Both the two-dimensional commutative 
and Moyal planes can be viewed as the limits of a fuzzy sphere of infinite radius. 

2. The two dimensional Model and its fuzzy version 

We are interested in the model: 

= Tr [a^'^ [L,-, [L;,0]] +Zp<I>2 + c<I>4] , (2.1) 

where <I> is a Hermitian matrix of size N, b and c we. mass and coupling parameters respectively. 
Li is the angular- momentum generator in the N dimensional unitary iiTcducible representation of 
SU (2). Since a rescaling of <I> will allow us to set a = 1, the entire phase diagram can be explored 
by ranging through all real values of b and positive values of c. The conventions of [Q] are a = 
b = arF?-, c = aXR^. 

The infinite matrix limit of the action can be taken and corresponds to a real scalar field on 
a round sphere of radius R and Euclidean action 

S[0] = j^^ d'-n {<p^^<p + r/?202 + XR^ip^) (2.2) 

where .J^^ = Ei=i,3-^-^ the usual angular momentum generators. 

The eigenvectors of [L;,[L/,-]] in (pTI) are the polarization tensors %„ (normalised so that 



^Tr(Fjjjy/,„) = 1) and it has eigenvalues /(/ + 1) with degeneracy 21 + 1. This is precisely the 
spectrum of the Laplacian on the commutative sphere truncated at angular- momentum A'^— 1. 

This particular model was chosen because of its simplicity. The diagrammatic expansion of the 
model ( |2.2[ ) has only one divergent diagram, the tadpole diagram, is Borel resumable, and defines 
the field theory entirely. In the fuzzy version, the tadpole splits into planar and non-planar tadpoles, 
which ai-e also the only diagrams that diverge in the infinite A'^ limit. Their difference is finite and 
nonlocal and is responsible for the UV/IR mixing phenomena of the disordered phase [^]. 

Even though the scalar field on either commutative or fuzzy spheres cannot have a phase 
transition, since they have finite volume or a finite number of degrees of freedom, phase transitions 
may be found when the matrix dimension or the radius of the sphere become infinite. 
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The fuzzy sphere can be recognized by introducing coordinates (Xi,X2,X3) proportional to the 
angular- momentum operators 



2R 



Li. 



They must satisfy the algebra 



Xf + Xi+X^=R% 



r 1 ® 
[Xi,Xj\ = i—EjjkXk 

K 



where 



2R^ 



is the parameter of non-commutativity, R is the radius of the sphere and 1 is the 



unit operator. The non-commutativity parameter depends on the matrix size, A'^, and the radius of 
the sphere, R. By taking different limits we can access different spaces: 
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2.1 Order parameters 

A suitable set of order parameters can be identified from the coefficients of a mode decompo- 
sition in terms of the polarization tensors basis |§] 



12 



(2.3) 



We have separated explicitly / = and 1=1 from the expansion to identify two observables whose 
expectation values we used to identify the respective phases. The observables are |rr(<I>)| and 

3 

p2 := p^,p^, = ^ {Tr{La^)f. The total power in all coefficients is given by Tr{^^) = ^{Tr{^)f + 

a=\ 

N(N^-i) P^ + k'mP ^^'^ used to estimate the importance of the neglected higher 

1=2 m=-l 

modes. 



2.2 The phases 

This model (2J_) presents three phases. As a generic illustration of their properties. Fig. [T](a) 
and Fig.|](b) show the dependence on the mass parameter b of the probability distributions of Tr{^) 
and p, respectively, for {a = 1, c = 40, A'^ = 4}. 

Disordered: Found for \b\ "small", the configurations fluctuate close to <I> = 0. This is con- 
firmed on the figure, < |rr(<I>)| >~< p >~ 0, but also < rr(<5>^) >~ (not shown). 

Non-uniform ordered: As \b\ increases, the figure shows multiple symmetric peaks for the 
probability distribution of Tr{(^) whose height decreases with increasing |rr(<I>)|, and multiple 



peaks not centered neai^ zero for p. Furthermore rr(<I>^) is much larger than both 



<|rr(o)|>- 



and 



< p > so that higher modes actually dominate. In particular, the most probable configuration is not 
rotationally invariant and we have spontaneous breakdown of rotational invariance. The probability 
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distribution of Tr{^) has A'^+ 1 symmetric peaks located approximately at — Ik)^ —h jlc where 
^ = 0, 1, . . . ,A'^, while the probability distributions of p and 5[<I>] have (A'^ + l)/2 peaks for A'^ odd 
and N /2 + 1 for N even. 

Uniform ordered: As \b\ becomes large, the figure shows two symmetric peaks for the prob- 
ability distribution of Tr{^) corresponding to the outer peaks of the non-uniform ordered phase 
and located approximately at Tr{^) ±Ny^—b/2c, but just one peak near zero for p. Further- 
more, < rr(O^) >~< Tr{^) >^ /N indicating that the power in higher modes is negligible. This 
is generic and indicates that sj —hjlc 1 and the rotational symmetry is thus restored. 



N=4 , a= 1 , b= -75.0, c= 40 N=4 , a= 1 , b= -75.0, c= 40 




(a) Probability distribution of Tr{^) (b) Probability distribution of p . 



Figure 1: Figures (a) and (b) show the typical behavior of the observables rr(<I>) and p in a region of the 
phase diagram where decreasing h passes the system through the three phases. 



3. Simulation and Results: The specific heat and the phase diagram 



We are interested in the phase diagram of the model ([2. Ij). To identify it, we used the coordi- 
nates in parameter space of the peaks of the specitic heat, C :=< > — <S >^. The relevant set 
of parameters is \N ^h^c] where h and c depend implicitly in R. It is possible to further reduce by 
one the number of parameters by finding a scaling {^,c} — > {hN^'' ^c'N^'^^ . If we find such and 
Q^, the model becomes independent of A'^ and automatically yields an infinite matrix limit. 

The simulations show that in the non-uniform ordered phase, the fuzzy kinetic term (propor- 



tional to a in (p. 1|)), is negligible compai^ed to the potential term (the other terms). There exists 
an exact solution for the corresponding limit of a = in the large A'^ limit called the pure potential 
model [0]. This model predicts a third order phase transition between a disordered and non-uniform 
ordered phase at c = /AN. Figure ^confirms numerically the convergence of the disordered/non- 
uniform ordered transition towai^ds this exact critical line of the pure potential model. 

Numerically, it is not difficult to find the coexistence curve between the uniform ordered and 
disordered phases which exist for low values of c. On the other hand, the coexistence curve between 
the two ordered phases is difficult to evaluate because it involves a jump in the field configuration 
and tunnelling over a wide potential barrier. 
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Figure 2: Plot of the specific heat at the disordered/non-uniform ordered transition for increasing and its 
N ^ °° limit, the exact pure potential model. 
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Figure 3: Phase diagram obtained from Monte-Carlo simulations of the model (2.1) 



The phase diagram obtained by Monte Carlo simulations with the Metropolis algorithm is 
shown on figure |[ That plot shows the phase diagram for the model ( |2.1| ). The data have been 
collapsed using the scaling form defined above with db = —3/2 and 6^ = — 2. It is remarkable that 
this scaling also works for the exact solution of the pure model potential. 
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4. Conclusions 

The numerical study showed three different phases. In one of those phases, which does not 
exist in the commutative planar theory, the rotational symmetry is spontaneously broken. The 
other two phases have quaUtatively the same character as the phases of this latter model. The three 
coexistence curves intersect at a triple point given by 

{b^,c^) = {-0.\5N^I^,0m^). (4.1) 

Those three curves and the triple point collapse using the same scaling function of A'^ and thus 
give a consistent N ^ °° limit. Thus, all three phases, and in particular the new uniform disordered 
phase, and the triple point survive in the limit. We will discuss these issues more completely in 
Ref. [f. 
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